FINAL REPORT ON 


Development of a Global Pollution Model For 
CO, CH^, and CHgO 



1983 East 24th Street 
Cleveland, Ohio 44115 


NASA Research Grant No. NSG 1030 
Grant Period: May 1, 1974 - August 14, 1974 


(NASA-CR-139484) DEVELOPMENT OF A GLOBAL N74-31115, 

POLLUTION MODEL FOR CO, CH4, AND CH20 
Final Report, 1 May - 14 Aug- 1974 

Cleveland State Univ.) 16 p HC $4.00 Unclas 

CSCL 04A G3/20 46362 





I 


Development of a Global Pollution Model for 
CO, CH^, and CH^Q 

K by 

L. K. Peters 

ABSTRACT 

The remote sensing of atmospheric pollutants from satellites will be 
feasible within the next few years. To effectively interpret these data, 
global pollution models that describe the transport and the physical and 
chemical processes occurring between sources and sinks of pollution are 
essential. This report describes the current status of a global pollution 
model for carbon monoxide, methane, and formaldehyde. 

The study considers the physico-chemical action of these three pollutants 
in the troposphere. This geographic restriction is convenient since the tropo- 
pause provides a natural boundary across which little transport occurs. The 
data on sources and sinks for these pollutants is based on available informa- 
tion and assumptions relative to the major man-made and natural contributions. 

The distributions and concentrations of methane, formaldehyde, and carbon monoxide 
in the atmosphere are interrelated by the chemical reactions in which they par- 
ticipate. A chemical kinetic model based on the pseudo-steady state approxima- 
tion for the intermediate species has been developed to account for these reactions. 

The numerical procedure being used to mathematically describe the pollution, 
transport is a mass conservative scheme employing an integral flux approach. 

It is fourth-order accurate in space which is desirable in simulating convective 
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processes in three space dimensions. Since computer storage places restric- 
tions on the scale of transport processes that are explicitly calculated, 
smaller sca]e mixing is described using an artificial diffusivity. 

At the present, the computer model is successfully functioning for short 
time integrations. Some additional work remains to develop a plotting routine 
and to generate a realistic modeling parameter data tape. 
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INTRODUCTION 

In the MAPS (Measurement of Air Pollution from Satellites) experiment 

** ( i \ 

proposed for:, NIMBUS 0 by.Reichle' the global distributions of carbon 
monoxide, sulfur dioxide, nitrogen dioxide, methane, ammonia, and aerosols 
would be measured over a two year period. The utility of satellite remote 
sensing data of pollutants .will be greatly enhanced if meaningful global 
models of pollution are developed. These models will be required to realis- 
tically analyze the atmospheric transport of pollution from its sources to 
sinks. Furthermore, these models can be significant in placing anthropogenic 
sources in proper perspective on a global scale. This report briefly describes 
the development and current status of a tropospheric global model for carbon 
monoxide (CO), methane, (CH^), and formaldehyde (Ch^O). 


The analysis is accomplished by geographically distributing the sources 
and sinks of CO, CH^, and CH^O, and simulating their convective and diffusive 
transport by a numerical solution on the computer of the turbulent diffusion 
equation. The atmospheric phenomena of these three species are coupled through 
atmospheric chemical reactions that occur. Thus, the three species must be 
considered simultaneously. Generally speaking, the oxidation of methane pro- 
duces formaldehyde which decomposes to carbon monoxide. Other sources and 
sinks of these pollutants are also operating. 


Other analyses to the present have utilized a global residence time 
(2 3 4 5 6 e a ) 

approach' * 9 9 5 * « Models incorporating a multiplicity of sources and 

sinks have not, to the author's knowledge, been attempted. The model outlined 
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in this study employs known source and sink strength data, the atmospheric 
chemistry of the pollutants in question, monthly averaged climatological data, 
and the turbulent diffusion equation for each species to establish global 
concentration distributions. ■ - > 

GENERAL MODEL DESCRIPTION 

The model development is restricted to the troposphere. This is a logical 
boundary since the tropopause provides a natural surface thru which the rate 
of mass transfer is relatively low. Furthermore, the photolytic decomposition 
of CO 2 appears to be unimportant as a source of CO in the troposphere^, and 
this enables one to decouple the CO transport from the CO 2 transport. 

Sources and Sinks of CH^ 

The sources and sinks of methane appear to be reasonably well understood 
at the present. The anthropogenic sources are largely the result of internal 
combustion engines and oil drilling and refinery operations. These emissions 
can be fairly well mapped based on automobile density and industrial activities. 

The natural sources apparently far out strip the man made sources - the 
principal ones being decaying vegetation and other biological action. Some 
of this biological action occurs within marine environments, and as a result 
the surface waters of the oceans, bays, and rivers appear to be supersaturated 
with methane. Lamontagne, Swinnerton, Linnenbom, and Smith ^ have reported 
equivalent surface ocean and sea water concentrations about 1.2 to 1.7 times 
the corresponding atmospheric concentrations. Specifically in open tropical 
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ocean waters, the surface concentrations (4.7 x 10 ml/1) corresponded to an 

equilibrium atmospheric concentration of 1.80 ppm v/hereas the measured atmos- 
pheric concentrations averaged 1.38 ppm. Bay and river waters appear to be 
even more heavily supersaturated. Their results specifically cited the following 
supersaturation ratios: Chesapeake Bay -14.3, York River - 21.2, Mississippi 

River - 5.67, Potomac River - 36.0. -These values may also be- affected by 
local pollution problems. The data of Brooks and Sackett^ on the coastal 
waters of the Gulf of Mexico generally support Lamontagne et al's results. 
However, they report that in the Yucutan area, where there is a major upwelling 
of deep water with low hydrocarbon concentration, the 'Gulf of Mexico acts as a 
sink for methane. 

The principal sink mechanism for methane appears to be in the homogeneous 
gas phase reaction of methane with hydroxyl radicals. 

CH 4 + OH ■+ CH 3 + H 2 0 (1) 

The methyl radical can subsequently undergo reactions which results in for- 
maldehyde and ultimately in CO formation. Thus, this sink for methane provides 
one natural source for formaldehyde and carbon monoxide. The following sequence 
of reactions is responsible for producing the hydroxyl radical 


0 3 + hv 0( ] D) + 0 2 (2) 

0( 1 D) + H 2 0 ■+ 20H . (3) 

0( ] D) + M 0 + M (4) 

0 + 0 2 + H 0, + M (5) 


It will be noted later that the hydroxyl and atomic oxygen (0) are also impor- 
tant in reactions with CO. 
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Sources and Sinks of C^O 

The anthropogenic sources of formaldehyde appear to be relatively small - 
the main onks being direct emission from automobile exhaust and formation 
during photochemical smog episodes. Jbe.se estimates cm be fairly reliably 
based on past auto exhaust emission estimates and .studies . _ 

The only apparent natural source for CH^O is from the methane oxidation 
just cited. Levy^’^ and McConnell, McElroy, and Wof sy^ ^ have suggested 
the following steps in the formation of formaldehyde by this mechanism. 


ch 4 + 

OH + CH 3 + H 2 0 

(1) 

ch 3 + 

0 2 + M + CH 3 0 2 + M 

(6) 

CH 3 ° 2 

+ NO -+■ CH 3 0 + N0 2 

(7) 

ch 3 0 2 

+ CH 3 0 2 ^ 2CH 3 0 + 0 2 

(8) 

CHjO •( 

■ 0 2 CH 2 0 + H0 2 

(9) 


Thus, the formaldehyde formation and concentration is directly coupled to the 
methane distribution. 

The main sink of formaldehyde is in photochemical decomposition and the 

i 

reaction with hydroxyl radicals. The following reactions appear to be impor- 
tanJ 10 ' 11 * 12 ). 

CH 2 0 + h V + CHO + H (10) 

CH 2 0 + hv -*■ H 2 + CO (11) 

CH 2 0 + OH 


CHO + H 2 0 


(12) 
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The production of CHO also leads to carbon monoxide formation via^’^ 

CHO + 0 2 + CO + H0 2 (13). 

Therefore, the source and sipk distribution of formaldehyde is primarily due 
to homogeneous gas phase reactions and is coupled to the methane and carbon 
monoxide distributions. 

Sources and Sinks of CO 

The sources, sinks, and concentrations of carbon monoxide have been 

(31 

summarized by Bortner, Kummler, and Jaffe v . The world-wide anthropogenic 
sources are estimated to be slightly in excess of 300 million tons/year with 
nearly two-thirds resulting from motor vehicle emissions. The remainder is 
distributed between stationary combustion sources, industrial' processing, and 
incineration. Therefore, these sources would be distributed largely according 
to motor vehicle density. 

The major natural sources of carbon monoxide appear to be the oceans, 

forest fires, terpene photochemistry, and gas phase reactions. The studies of 

(5l (151 

Junge, Seiler, and Warneck' , Seiler and Junge' , Swinnerton, Linnenbom, and 
Check^^, and Lamontagne, Swinnerton, and Linnenbom^^ indicate that the 
level of excess CO in the ocean corresponds to an equilibrium air phase concen- 
tration of about 3.5 ppm. These data were obtained during ocean cruises. It 
is, of course, reasonable to expect that the river, lake, or ocean regions near 
urban areas, where the CO concentrations may be considerably larger than 3.5 ppm, 
may act as a sink for CO. Furthermore, it is possible that the oceans at the 
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high latitudes serye as a .sink for CO produced by oceans at the low latitudes > 
since the warmer tropical waters would likely have a higher biological activity 
producing mo,re CO. With the warmer water, the solubility of CO is reduced. 
Transport over colder waters with greater CO solubility capacity then creates 
the possibility of these sections acting as a sink for CO. Thus, it is plausible 
that the oceans act as both a source and sink for CO. Similar arguments could 
be made for methane. 

In addition to terpene photochemistry and forest fires (combined sources 

are estimated at 23 x 10^ ton/year^^ ) , the other principal natural source of 

CO appears to be the gas phase reactions cited earlier. However, these cannot 

be divorced from the reactions which consume CO, of which the following seem to 
( 3 ) 

be important' ' . 


CO 

+ 

OH -> C0 2 + H 

(14) 

CO 

+ 

0 + M C0 2 + M 

(15) 

CO 

+ 

N 2 0 Sur l ace co 2 + n 2 

(16) 

CO 

+ 

H0 2 + C0 2 + OH 

07) 


Reaction 21 ,22,23) repor tedly first order in CO but zeroth order 

in N^O. This is a surface catalyzed reaction and requires, for complete 

accuracy, detailed information on the atmospheric aerosol as to size distribution 

and chemical composition. Reaction (17) has been suggested by Westenberg^) 

to be important in atmospheric pollution problems. However, the results of 

(25) 

Davis, Wong, Payne, and $tief v { indicate that it is unimportant in the overall 
oxidation processes of CO. As a result, the current model ignores Reaction (17). 
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Another natural sink of CO of seemingly large significance is the soil. -* 

The recent work of Inman^^, Inman, Ingersoll, and Levy^^, and Ingersoll 
( 28) 

and Inman' 1 point up this significance. Recent field studies in which soils 

^ i 

-3 

were -'exposed, in situ to test atmospheres containing 115 mg m” CO (100 ppm)*.' t • 

- 1-2 

showed average uptake rates that varied from 3.5 ug s m for desert areas to 
-1 -2 

16.5 yg s m for tropical deciduous forest areas. By using an average con- 

centration driving force of 57 mg m CO (50 ppm) and assuming conditions of 

atmospheric pressure and 293 K, the mass transfer coefficient at the surface 

corresponds to 0.594 x 10~^ m s""* and 2.38 x 10"^ m s""*, respectively. These 

values are probably indicative of the rate of the biological reaction that is 

occurring near the surface of the soil and may simply be representative of the 

type and concentration of the soil microorganisms utilizing CO. One of the big 

uncertainties relative to soil scavenging is the determination of which fraction 

( 29 ) 

of these microorganisms are anaerobic methane-producing s ' and what fraction 

(15) 

are aerobic CO^ producing. Furthermore, Seiler and Junge v ' have suggested 

-3 

that at low concentrations (around 0.23 mg nT at 298 K) a temperature depen- 
dent equilibrium of CO. above soils occurs. If this is truly the case, then 
the soils can act as either sinks or sources for carbon monoxide in much the 

i 

same manner as do the oceans. 

\ 

A source and/or sink common to all three species is leakage from and/or 
to the stratosphere. Since the initial model is restricted to the troposphere, 
this leakage is considered as sources and/or sinks. In the specific case of CO, 
leakage from the troposphere to the stratosphere is reasonable to expect since 
the CO that escapes thru the tropopause will typically undergo chemical reactions 
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and not return to the troposphere as CO. This is substantiated by the vertical 
■profiles of CO which show a. decrease of CO mixing ratio .with height above the 
tropopause. 

The leakage of CH^ thru the tropopause would be similar to that for CO 

since vertical profiles of CH^ also show a decrease with increasing altitude 

f 13 311 

above the tropopause^ 5 '. This type of information has not as yet been 

found for formaldehyde. 

CURRENT STATE OF MODEL DEVELOPMENT 

The general approach in simulation of the CH^ - CH 2 0 - CO cycle consists 
of the following steps: 

a. Initialize the CH^, CH^O and CO concentrations in the troposphere. 

In order to conserve computer time, the initial concentrations 
selected are approximately that expected in the atmosphere. This 
should not affect the final results but simply the time required 
to reach steady state. 

b. Distribute the sources and sinks of the various species on the 

t 

Earth's surface and at the tropopause consistent with the 
physico-chemical considerations. This involves proper interpre- 
tation of oceans and lands as sources and/or sinks of the par- 
ticular species. As a first approximation, the tropopause is 
being considered as a zero flux boundary. Thus, all pollutant 
generation and consumption is entirely within the troposphere. 
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c. Solve the three coupled unsteady state turbulent diffusion 
equations for CH^, CH^O, and CO with the boundary conditions 
established by b. The appropriate climatological data^^ 
is used to establish the wind field, temperature field, and 

water vapor field. The coupling of the diffusion equations results 
from the gas phase reactions creating homogeneous generation terms. 
This, therefore, accounts for the chemical sources and sinks present. 
A simplified chemical kinetic model based on the pseudo-steady 
state approximation for the intermediate species has been developed 
to account for these reactions. 

d. Continue the integration in time until reasonably steady concen- 
tration distributions exist. 

The inherent advantage to using this procedure is that one does not 
presuppose the atmospheric concentrations of the three pollutant species 
being studied. This, therefore, provides a meaningful test of the distri- 
bution of sources and sinks. This should not imply however, that there are 
not uncertainties present. But these uncertainties are mostly associated 
with the strengths of the sources and sinks, and a primary goal of satellite 

i 

remote sensing is to establish the magnitude of these source and sink strengths 

At the present time, the computer program to describe the physico-chemical 
behavior of the three species is successfully functioning on the computer 
system at NASA Langley Research Center for short time integrations using res- 
tricted modeling parameters. Since the space length scales in the numerical 
calculation are quite large, the simulation must be continued for relatively 



12 


long periods of time. Current results indicate that real time simulation on 
the order of ten to thirty- days is required. Furthermore, computer time to 
simulation time appears to be on the order of 1:10 to 1:20. Therefore, simu- 
lation with one set of modeling parameters will likely require twenty to thirty 
hours- of CDC6600 time. As a result, the possibility of using the NCAR CDC7600 , 
which is approximately five times faster, is being investigated. 

To supplement the main calculating program, a data tape incorporating 
realistic modeling parameters and a computer routine for plotting isopleths of 
the calculated concentrations are being developed. The successful coordination 
of the three separate phases should be established before any long term numeri-~ 
cal integration is attempted. It is hoped that these studies can be con- 
tinued. 
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